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ABSTRACT 

Numerical integration plays an important role in satellite orbit 
determination. This paper presents the general philosophy of numeri- 
cal integration, a description of the often used multistep numerical 
integration algorithms pertinent to orbit determination, and the deri- 
vation of the formulas and their various forms used in these multistep 
algorithms. The coefficients for different forms of these formulas are 
presented in rational form up to order fifteen in the appendix. 
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COWELL TYPE NUMERICAL INTEGRATION 
AS APPLIED TO SATELLITE ORBIT COMPUTATION 


by 

Jesse L. Maury, Jr., and Gail P. Segal 

Goddard Space Flight Center 


INTRODUCTION: GENERAL PHILOSOPHY 


Many problems involving ordinary differential equations cannot be solved explicitly or 
analytically. It is for this reason that numerical techniques for approximating solutions of such 
equations were developed. The advent of high speed computers which can handle the tedious 
arithmetic involved has made these techniques even more attractive and useful. Using a compu- 
ter, it is possible to extend these numerical techniques to a degree of precision far higher than 
any hand calculation could ever achieve. 


Of particular interest are the discrete variable methods which yield approximate solutions of 
the problem y f (x, y) at a set of discrete points x, x + h, x + 2h, . . . where h is the step size. 

In general, the discrete variable methods applied to initial value problems can be classified as 
either one-step methods or multistep methods. The one-step methods require knowledge of the 
value of the function at only the previous point while the multistep methods require this knowledge 
at a certain number of preceding values. That is, to approximate the value of the function at x + h, 
a one-step method would need only knowledge of the value of the function at x while a multistep 
method would require this knowledge at the points x, x - h, x - 2h, x - 3h, . . ., x - nh. 

At first, one might think that the one-step methods would be more advantageous in obtaining 
the approximations since they require only one previous value, one backpoint. However, the error 
committed in using the formulas of any one-step process over a given interval is generally larger 
than the error incurred in a multistep method. Also, to go one step forward with a one-step 
method requires more evaluations of the function, and, in the multistep method, increasing the 
order (the number of backpoints used) does not necessarily require a concomitant increase in 
evaluations. Furthermore, since large orders of a multistep method are easily attained, multistep 
methods are highly accurate with relatively large increments of the independent variable. 

In the realm of orbital dynamics, the use of numerical techniques is virtually dictated. It is 
almost impossible to solve analytically (i.e., explicitly) those equations which represent the mo- 
tion of a satellite. Analytical solutions such as Brouwer or Two Body Motion are sometimes 


1 





employed, but at best they use only limited approximations of the real forces which affect a 
satellite’s motion. With the numerical approach, the expression^ of these forces do not have to be 
truncated after the first few terms: they can be expressed in their entirety. 

Some of the computer programs which use numerical methods to compute the motion of 
artificial satellites are: 

D.O.D.S. — Definitive Orbit Determination System 
May 15, 1968 

Space Systems Analysis and Computer Programming Services 
Contract NAS 5-10022 

Prepared by 

Scientific Satellite Systems Department 
Federal Systems Division 
International Business Machines Corporation 
Gaithersburg, Maryland 

Noname — An Orbit and Geodetic Parameter Estimation System 
Aug. 1968 

Contract Number NAS-5-9756-71D 
Prepared by 

Wolf Research & Development Corporation 
Applied Sciences Department 
College Park, Maryland 

Prepared for 

Mission and Trajectory Analysis Division 

National Aeronautics and Space Administration 
G.S.F.C., Greenbelt, Maryland 

Lungfish — Lunar Gravitational Field in Spherical Harmonics 
Feb. 1966 

Contract No. NASI -4998 

Prepared for the Space Mechanics Division of the Langley Research Center 
Prepared by Computer Usage Company, Inc. 

Trace — Trace -C Powered Flight Trajectory Determination Program 
May 1965 

Report No. TOR - 4 69(5352)-! 

Prepared by Aerospace Corp. — 

C. S. Christensen, A. R. Jacobsen and R. J. Mercer 

This paper describes how multistep numerical integration is started with a one-step process, 
exemplified by the Runge-Kutta method; how the multistep process is used in orbit determination, 
exemplified by Cowell type formulas; and derivation of predictor and corrector formulas for 
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equations of the first and second orders. Also included, in the appendix, are the coefficients for 
the multistep methods discussed in the text. 

In the discussion, y and f are 3 -space vectors. The independent variable is x, while 

I y I = (yi 2 + y 2 2 T y 3 2 )^ • 

DESCRIPTION OF INTEGRATION METHOD 
I Starting the Multi step 

The multistep numerical integration method of solving differential equations requires a 
knowledge of preceding values (backpoints). Consider the initial value problem 

y 1 = f (x. y (x)) 

y(*o) = y 0 * 

We need to know the values y (x^ = y v y (x 2 ) = y 2 , . . y (x^) = y mM ,y (x^) = y m where - x * h, 
x 2 r x + 2h, . . x m ^ 1 ~ x + (m - 1) h, x m = x + mh, h being the step size. These va i ?* s are needed 
to determine from evaluation of y' - f (x, y(x)) — more simply written f (x, y) — the backpoints y^ , 
yj- . . • m yi» y 1 , y 0 required by the multistep algorithms. (In physical terms, this may be con- 

sidered as having for each x x a position y. and a velocity .) 

To produce the initial backpoints used to start the multistep process, a one-step numerical 
integration method such as Euler’s method, Taylor T s expansion, Runge-Kutta, etc., is used. Each of 
these methods requires a knowledge of only one preceding value of y(x). Thus the initial value 
y(x 0 ) = y 0 is sufficient to initiate the one -step 11 starter” for a multistep process. 

A commonly used one-step method is the Runge-Kutta which computes y t , y 2 , . . . as follows: 
Given the initial value problem 

y 1 = f(x,y) 

y( x o> = y Q * 


The formula used is 


y nM - y n + J( k i + 2k 2 + 3k 3 + 4k 4 ) 

n = 0.1,2,.-. 


where 


k l = hf ( X n’V n ) 
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/ h k i' 

k :i = hf \*„ + T' y n ' T. 


k . - h f ( x + h , y t k , ) 

As can be seen from the above equations, a fourth order Runge-Kutta process requires four 
evaluations of the derivative y’ = f(x.y > for each step forward. 

By way of remark, the following should be considered. Applying this Runge-Kutta process to 
each of the three (usually complex) equations of motion of a satellite to produce position and 
velocity coordinates is inefficient. Furthermore, to achieve the required accuracy necessary in 
orbit determination analysis, the step size h must be very small. The error incurred by this 
fourth order Runge-Kutta is of the order h 5 while the corresponding local error for a multistep 
process is of the order h PTl where P is the order of the multistep method which is usually higher 
than 4. Thus, the step size of the Runge-Kutta starterums^ be a fraction of the step size of the 
multistep process. This is an important consideration in programming the multistep algorithms. 

There do exist multistep methods used as starters. These methods employ a time-consuming, 
iterative procedure to produce each backpoint and it is questionable whether they are more effi- 
cient than the one-step methods. In any event, the time required to set up the starting table of 
initial backpoints for the multistep process is usually a fraction of the total computation time. 

Any gains in efficiency accrued by these iterative schemes are, at most, marginal while the sim- 
plicity of the one-step methods make them desirable. 

II. The Multistep Algorithm 

Assuming now that for the initial value problem 

y’ = f(x. y) 

y(* 0 ) r y 0 


= hf *„ r ~r- y 


we have generated the backpoints yj; , y^.j, ... yj, y ‘ , y<J by some single step process. (We may 
write y‘ = f 0 , y ; = i v ..., y^ = f m to mean y ' = f(x 0 + ih. y(x 0 ; ih)).) With this set of back- 
points, y 0 , y 1( y 2 , . . . y , y m , the multistep process can be started. These values are used in 
an extrapolator or predictor to compute y m + I . The predictor considered here is the Adams- 
Bashforth (Henrici) which has the form 






4 



where represents a difference operator (dis- 
cussed later) operating on and employing the 
backpoints y m , C,, .... y'. nt2 , y^ M . 

The predicted value of y mM is used with 
x n+] to evaluate 

y' - y) 

for y„* i - Thts value of y^ M is then employed 
in a corrector formula which yields a new 
value for y m+ : . The corrector discussed here 
is the Adams-Moulton (Henrici) which has the 
form 




y. T h {“o '°y: + i T ar^y^i 


+ a. 


2 y ' * . T ... + a * V n y 1 ) 

m l n ■ 7 m t 1 / 


We now have 


two 

Pv 


values for 
and 


y . . . d 

predicted value, say p y r>+1 , and a corrected 
value, say c y m + 1 . These two values are com- 
pared. If the absolute value of their difference, 

I c y tn + 1 “ p y m + 1 S > is not less than a given toler- 
ance, the Cl y m+J is used (i.e., substituted for 
p y f11 + j ) with x m+I to again evaluate f (x. y) for a 
new value of yJ| M . The corrector is then used 
again with this new value of y' m + l to calculate a 


new l 
rector i 

where c * 


This iteration process on the cor- 


repeated until 

- p 


+ i 


y t „ 


- < 


1 y n + l ' 5 

meets the tolerance. A 




simple flow chart may describe this more clearly. See Figure 1. 


Figure 1 — Predictor-corrector algorithm applied to the 
initial value problem y* = f(x,y), y(x Q ) = y Q . 


When the iteration process has converged (i.e., the criterion on | Cl + l y m+1 - C1 y m + l | has been 
satisfied), the final computed value for is entered in the backpoint table. Then, where the 
points y ’, y ;, ..., y^j, yj; were used to determine y^ +1 , the points y J, y^, . . y ;, y n + 1 are now 
used to determine y' m + 2 . Etc. 

Note that in the Adams-Bashforth predictor, no knowledge of the value y m+1 being derived is 
needed while such knowledge (namely a value for + I ) is needed in the Adams-Moulton corrector. 
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Equations like Adams -Moulton corrector (closed form equations) have smaller truncation errors 
as well as desirable stabilizing characteristics. The predictor is used to obtain an estimated 
value for y m+1 good enough to keep the number of corrector iterations low. This predictor - 
corrector algorithm is well known and it has been shown by various authors that for a sufficiently 
smaD step size, h, the successive corrected values obtained converge to the unique solution of the 
closed form equation provided the function being numerically integrated is sufficiently smooth. 

The above discussion considered numerical calculations for deriving values of y (and con- 
comitantly y 1 ) at discrete points from the initial value problem 

y* - f ( x. v ) 
y ( x 0 ) " y 0 • 

The same technique could be used on any initial value problem of the form 


y' r '> = f (x, 

y (n ~ u (x 0 > - y 0 (n ' u 

to solve for y (n_,) (x i ). In particular, we are interested in calculating y j +I from the backpoints 
y”, y^'_ , . . . since, in general, satellite orbit determination involves the initial value problem 

y" = f (*. y. y' ) 

y'( x 0 ) = yo 
y( x 0 ) = y 0 

This could be approached by generating an initial set of backpoints for y" and y then using 
y”, y”.,, ... to calculate y^, and using y„; , y^.j , ... to calculate y m + 1 employing the same tech- 
nique described above in both steps. However, certain advantages accrue if we use a mathemati- 
cally equivalent technique which derives y m< 1 directly from the backpoints y ”, i ■ • • ■ For 
one, it is necessary to keep only one set of backpoints — the retention of y^, y^_j , ... is obvi- 
ated. Secondly, we often must work with the problem 

y" = f(x.y) 


y( x 0 ) = y 0 
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when only conservative forces are involved (i.e., no drag or other energy dissipating forces). In 
this situation, when y “ +1 has been satisfactorily determined, y^ M can be calculated by evalu- 
ating the corrector 




* T 70 t 


a.-V. 


L ‘ V" 


Ci} 


only once. 

Consider now, working with the initial value problem 

y" " f (x, y) 
y’ (x 0 ) = y ' 0 

y( x o> = y 0 


Here, the predictor -corrector approach is the same. The difference exists in the polynomials: in 
particular, the coefficients are different. The formulas considered here are generally referred 
to as Cowell type formulas. They are: 


Stormer Predictor 


y.*i 1 - y-t + h2 + /W + 


Cowell Corrector 


y = 2y - y , + h 2 ( >3 * y 1 ' +/3*V i y , ‘ + /3*V 2 v " + +/3*V n v M \ 

y m+l y m-l 1^0 y m+l ' 1 y m+l H 2 y m+] y ot + 1 J 


In the most general form of the initial value problem 

V M = f (x, y, y' ) 

y‘ (* 0 ) = yo 
y(* 0 ) = y 0 - 

yJ, M is derived from the backpoints y", y^, , . . . using the Adams formulas while y m+1 is 
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obtained from the same backpoint set using 
the Cowell formulas. In testing for con- 
vergence of the corrector formulas, the sum 

! c, m * - c . y ' i + ' Ci + : y - c i y is com- 

pared to the tolerance. A flow chart of the 
process is given in Figure 2, 

HI. Derivation of Multistep Formulas 

These foregoing techniques are referred 
to as numerical integration. This appellation 
originates from the derivation of the methods. 
Consider again 

y' - f(*. y) 

y(* 0 ) = y 0 • 

Integrating both sides between and x n+! , we 
have 





or 



C x m + ; 



' y m r 

y (s)ds 

STORE TIN, 

aL value of 

J 

x m 

y^ +) N BACKPOINT TABLE 

v.m = y„ + J 

f *m+ 1 

| f ( s ) d S 

X TTl 


3 


where f(s) denotes f (s. y(s))- 

By replacing f (s) by a Newtonian type interpo- 
lating polynomial and integrating, it is possible 


Figure 2 — Predictor-corrector algorithm applied to the 
Initial value problem y" = f (x, y f y * ), y 1 (x 0 ) = y' Q , 

y (* 0 ) = Yo- 


to derive the Adams type polynomials which are used to approximate the expression 



f ( s ) d s . 


The error generated by replacing the function being integrated with a polynomial which is, 
effectively, integrated is usually obtained by integrating the local error associated with the interpo- 
lating polynomial. For example, it can be shown (Henrici) that the local error expression for 
formulas of the above type is of the form 


R r - C h pM y (p+1 ' (e) 
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where p is the order of the method, h the step size, f is a value between the largest and smallest 
values of x on the interval (x p . x p+1 ) , and c is a constant specific to the formula. 

The Cowell type formulas can be derived by a double integration of y" = f (x,y) and again 
employ’ng a Newtonian type interpolating polynomial (Henrici). These derivations are complex. 

A simp±er approach using difference operators avoids much of the difficulty involved in integrating 
the interpolating polynomials. This is the derivation given here. Using this approach, the opera- 
tor definitions lead naturally to the Adams -Moulton corrector. It is derived first. The other 
formulas follow easily from this derivation: first, the Adams -Bashforth predictor, then the Cowell 
corrector, and finally the Stormer predictor. 

In the ensuing derivations, some confusion may arise between the subscripts m and n +j. The 
predictors are derived for y mM , the correctors for y^. This is of no real importance since the 
same backpoints can be labelled either as y , y v or v , v . v .... 


A. Preliminary Definitions and Relationships 

In order to derive the formulas for multistep numerical integration, it is useful to develop 
several tools. Consider the following difference tables (Figures 3 and 4). The first column is 
formed by defining the values t(x + ih), i - 0, l, 2, ... for forward differences and f(* - ih) 
for backward differences. The second columns are formed from differences of successive values 
of the first column. The third columns, from .differences of the second. And so forth. (In both 
tables, the subtrahend is the value above the minuend in each column.) 


'f(xL 


f(x + bP'f(x)^ 

M f(x - 2hT^2f(x + h) + f(x) 

f(x r 2h) - f(x + h) + 30- 3((x + 2h) - 3f(x + h) - f (x) 

f <* + 2h) f(x . 3h) - 2 f (x + 2h) + f(x + h) 

f (x - 3h) - f(x -*• 2h) 

f(x + 3h) 

Figure 3 — Forward difference table. 

f (x - 3h) 


f (* — 2h) — f (x — 3h) 


fix - 2h) 


f(x - h) - 2f{x - 2h) + f(x - 3h) 



(x) - 3f (x - h) + 3f(x - 2h) - 3f(x - 3h) 


Figure 4— Backward difference table. 
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From these tables, we derive the following operator definitions: 


Forward Difference Operator (delta) 

£ f ( x) = f ( x - h) - 1 : x) 

A 2 f ( x) = A (A ( ( x) ) = f (x + 2h) - 2 f f x - h) * f (x) 

n 

A r f ( x) = A (A"" 1 f ( x) ) = J~Vir (”) l(x - (n - i ) h) 


Backward Difference Operator (nabla) 

(x) ” f(x)-f(x-h) 

V 2 f(x) - V(?f( x >) - f(x) - 2f(x - h) - f(x - 2h) 


V" f ( X ) 


n 

v(yn-l f(x) ) = (") f(x-ih) 

1 = 0 


These definitions simplify our difference tables. See Figures 5 and 6. 


"f(x) 



Af(x + h) £ 3 f M 

f(x + 2h) 4 2 f(* + h) 

Af(x + 2h) 

f(x +■ 3h) 


f (x - 3h) 


Vf(x - 2h) 

f(x - 2h) ^ 2 f (x - h) 



^f(x) 


Figure 5 — Forward difference table written in Figure 6— Backward difference table written in 

forward difference operator notation. backward difference operator notation. 


In addition to the difference operators, we define: 


Identity Operator 



Shift Operator 


Ef(x) f(x-'h) (4) 

f (x) = f (X + tj h ) 

( r) may be any real number) 

Differential Operator 

Df(x) = f(x) (5) 

D n f(x) = f^(x). 

On these operators, we define an algebra where, for any two operators L k and l 2 , L l ± L 2 
means the results of l 2 operating on f ( x ) are to be added to or subtracted from the results of i x 
operating on f t x) ; while multiplication, L : times l 2 , means l x operating on the results of L 2 
operating on f(x). For example, 

I f (x) - E " 1 f (x) = f(x)-f(x-h) - Vf (x), 

AVf(x) = A[f (x) - f (x - h)] 


= M(x) - Af (x - h) 


= f(x + h) - f ( x) - [f (X + h - h) - f(x - h)] 

- f(x + h) - 2 f (x) •+ f(x - h). 

It can be shown (Hildebrand) that these operators follow the laws of commutivity, associability, 
and distribution. 

With these definitions, we derive the relationships 

V = I - E" 1 (6) 

e = (i - V)' 1 = (v 

A = E - I. (8) 


Then from Equations (7) and (8), 


A = E - 


I I ~ I 2 + I ^ 

I - V ~ 1 “ I - V 
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But, i 2 


I and iv = 7. Hence 


A = 


I - v 


( 9 ) 


In addition to the above operator definitions and relationships, we need the series representa- 
tions for e x 9 j ^ , i - x > and C 1 “ x), and formulas for series multiplication and series 
division: 


X x* 

e x 1 + — . t 


1 ! 2 ' 


Zfr 


( 10 ) 


1 - X 


X ‘ X‘ T X 


2 r X 3 ^ 


I «- 


( 11 ) 


1 ~ X 


1 f X + x 2 


£ 


(12) 


X 2 X 3 V" 

•log(l-x) = X +— . . . - X / — — ■ • (13) 


For series division and multiplication, let the series Sj and s 2 be the arguments of the 
operation and s 3 the result. We define 


1 + a.x + a„x 2 




where 


a o = 1 


s, “ 1 + bj x + b 2 x 2 + **■ 



where 


b 


o 


1 , 


and for the resultant series s 3 - s 4 s 2 or :> 3 r s t s 2 we desire 


= 1 
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where 


c o 1 


Then, 

Series multiplication is defined as 


s 3 r 1 * ( b i + a i) x + ( b 2 T a i b i T a 2 ) x2 + ( b 3 r t a 2 b ] r a 3 )> 


i> (£‘-A 

L = 0 \ > =0 / 


(14) 


where 


a o = b o " 1 


and, 


Series division is defined as 


S l /S 2 ' S J “ 1 ' ( a i " b l) X + [ a 2 “ ( b l C l ' b 2>] x2 * [ a 3 “ ( b l c 2 + b 2 C 2 " b 3 )] 


X J + 


- 1 + 




(15) 


where 


C 


0 


1 . 


Note that series division is a recursive definition requiring c 0 , c,, c 2 , . . c n-1 to compute the n th 
coefficient, c n , of the n th term of the s 3 series. Note also, where Sj - 1, series division reduces to 


where 


since a. = o for i > 0 . 



1 


(16) 
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B. Derivation of Formulas 


\ 


Consider now the Taylor's expansion of an interpolating polynomial 


p(x + h) 


P(x) 


h 

T r 


p (,) (x) 



p< 2 ) (x) + 



P (r) (x). • 


Using the shift operator Ep(x) - P (x + h), the differential operator ir P (x) = P (n) (x), and the iden- 
tity operator l P (x) = p(x)> we have 


Ep(x) 


_h_ 

1 ! 


-D 2 ■* 



p(x). 


(Note that this is a finite expansion for any given n since P (x) is a polynomial, hence has only n 
derivatives.) 

Then, by Equation (10) the expansion of e x , 


E = e h D 


or, by relationship (7) is 


(I - V )- 1 = e hD • 


Taking the log of both sides, 


- log (I - V) = hD 


or 


hD 

- log(I " 


Multiplying both sides by v, 


57 = h [ -10,(1 - v> ] “ 

and employing Equation (13), the expansion of -iog( l * x), we have 
14 
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V = h 

n 

<i 

t 

O 

II 

X 

I 


y -2i_ 

/ _ ^ 1 + 1 
_ i = ° _ 

a 

Z— » i + 1 

_ i =° _ 


which by series division (16) is 


r? - 



(18) 


where n is the order of the interpolating polynomial and 


i; = l, 


• - f 1 

11 z_. j + 1 


i = i 


(19) 


This is the Adams-Moulton Corrector. Some of the coefficients, a.", are given in Table 1. 
For i = 0 to i = 15, see Table 2 in the appendix. 

Applying this to our initial value problem 

y n = f(x. y,y‘ ) 
y l ( x c> = y'o 

y(* 0 ) = y 0 


Table 1 


Coefficients of Adams -Moulton Corrector. 

i 

0 

1 

2 

3 

4 

5 

• 

a. 

i 

1 

1 

" 2 

1 

” 12 

1 

‘24 

19 
‘ 720 

, i 

3 

~ 160 


to obtain a corrected value, c y^ , when and y m have been predicted, an approximation of y ” cal- 
culated, and the other n - l backpoints y ^ 2 , . y^'_ nM determined, we have 


■ y m .i 


“ -}c 


15 




or 


y* = y«.i + h {^' " 1 [y» ■ y»'-i] 

12 f y ™ "" ~ y *>- 1 ’ ym-2] 


24 fy= " 3 C-1 + 3y m'-2 + y^-3] 


>•;: - + ( n 2 )y:., - Q y :-3 * 


(20) 


We now wish to develop the Adams -Bashforth predictor. Consider again Equation (17) and 
multiply both sides by relationship (7) noting that v& = Then 


T E - A = h 


( I - v )" 1 V 
“ log(I - V) 


D - h 


I - V 


log (I “ *) 


Now, employing Equations (11) and (13), the series representations respectively for and 

- log ( l - x) , we have 


A - h 


oc 

1=0 

n -- u 

i 

l 

^ 1 

t>\^ 

1 

u n 

*- | 

L.— 


which by Equation (15) series division is 


_ n 

L - 


V 1 


L i=o 


( 21 ) 


16 
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where n'is the number of backpoints (L.e., the order of the method) and 


1 -■•••» 


i = l 


(22) 


Some of the a are given in Table 2 . These coefficients are given rational form for i = o to i - 15 


in Table 2 of the appendix. 

Note that the derivation involved infinite 
series. However, since these operator rela- 
tionships are valid for polynomials, the corre- 
sponding series are finite. Hence, there exists 
n such that - 0 for all 1 > n. 


Table 2 


f Coefficients of Adams* 

-Bashforth Predictor. 


0 

1 

2 

3 

4 

. 

Cl 

1 

1 

2 

5 

12 

3 

8 

251 

720 

95 

288 


Thus, 








5 ^ 

~ 7^ r 4 V3 + 
12 8 



or 




- y m 


+ h 


y *' + 5 [y; 1 - y^'-J 


12 


[y« - 2 y^'-i + y-L a ] 


17 




~ + ^'-2 ' ^- 3 ] 


+ 



+ a. 

n 



(-D n y; 1 - 



(23) 


As previously noted, we have the problem of calculating from the backpoints y"_ lf y ^_ 2 , 
.... To achieve this, consider once again Equation (17). By squaring both sides we immediately 
have a formula involving D 2 y - y M . 


V2 


V - 

log(I - v) 


D 2 , 


(24) 


It is possible to obtain an expression for 


merely by squaring the series repre- 


sentation for 


V 


log (I “ V) 


log(I - V)J 

. However, a more suitable formulation can be derived as follows: 


Consider 


[-log(I - V)] 3 = D~‘ D [-log (I - '•)] 


where D' 1 is the informal integration operator (Hildebrand), the inverse of the differential opera- 
tor. Then 


D _I D [ — log ( I - V)] 2 


D -i 2 ['log (I ~ 


= D' 1 



from (13) 




from (11) 
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= D“ 


X 

L 


2H 


i + 1 




where 



k = 0 


Then, by integrating (i.e. using the operator D' 1 ), 


[-log(I - V)] 


L 


2H. V« t2 

i * i 

i + 2 



i '0 


(25) 


Using this expression in Equation (24), 


V 2 


h 2 


V 2 


2 V 2 


CD 

L 

, = n 


H ,+. 7i 

j + 2 


D 2 


which by series division (15) is 


V 2 


h 2 



D 2 


(26) 


where 
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and 



j = i 


( 27 ) 



W- 1 


Table 3 

i 

| Coefficients of Cowell Corrector 


1 

0 

1 

2 

3 

4 

5 

- » 

1 

! - 1 

1 1 

1 

0 ! 

1 

1 

i 

12 

" 240 

" 240 


and n is the order of the method. This is the 
Cowell corrector. Some of the coefficients, 
/5', are given in Table 3. For £*, i = 0 to 
i = 15 , in rational form see Table 4 of the 
appendix. 

Thus, 


V 2 y - y - 2 y . ~ y 0 

* sr ^ m - 1 ^ m — 2 






+ 0 


V7 3 _ 


1 

240 




or 




1 1 

m - 1 


6y " „ - 4y M 

^{ti - 2 - 


1 


y* ' (i) y«-i + (2) y -2 ' (3) y - 3 + + <' 1)n y ™’-n 


} 


(28) 


20 




As in the case of Equation (19) we need an extrapolator or predictor. This can be derived in 
the same manner as Equation (21), only this time, multiplying Equation (24) by relationship (7), 


V J E =■ h : 


NsM (rh) 


D 2 . 


Using Equations (25) and (12) 



which by series multiplication (14) is 


1 h2 Z] ^ ^ ^ 29) 

i e 0 


where 


^0 


1 


and 



2H 


j *i 


' Pi . 


(30) 


This is the Stormer predictor. Several of the 
coefficients, /5 i? are given in Table 4. For p. 
in rational form for i = o to i = 15 , see Table 3 
of the appendix. 

Thus, 


>m+ 1 


- 2y 


Table 4 


Coefficients of Stormer Predictor 


i 

0 

— 

1 

2 

3 

4 

— 

5 

3 

1 

0 

1 

1 

19 

3 1 

' 1 


12 

12 

240 

.. 

40 I 


= j I t 07 + + Ty3 
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or 



In recapitulation, we have derived the following formulas for numerically solving at discrete 
points the initial value problem 


y M = f (x, y. y ' ) 


y'( x 0 ) = 


y( x 0 ) = Vo- 


The Adams-Bashforth predictor 


= 


n+ 1 


n 



i =0 


where 


and 
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which is used to produce a first approximation of y^ + 1 for iteration in the Adams -Moulton cor 
rector 


n 



1 =0 


where 


= 1 


and 


i^O 


and the Stormer predictor 


= VJ E.v m = y.M - 2y m - y m _, = h 2 £ V y“ 

1 =0 

where 

4> = i 

and 



k s 1 


which produces a first approximation of y m + ] for iteration in the Cowell corrector 


: y. - ay.., + y m _ 2 = h 2 22 p; y” 

1=0 


where 


Q* 

H o 


1 
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and 


y-* 2H ,*. 

“ zL i + 2 

j = i 


,.6 *_ 


m 



k = 1 


C. The Summed Form 

It has been established (Henrici) that algebraic equivalents known as the summed forms of the 
foregoing equations considerably reduce the propagation of round-off error. These summed 
forms can be derived by defining the operators v -1 and v “ 2 as the inverses of T 1 and v* 

7-17 i l 7' 272 :: J 


and defining 


(32) 


y-i 


v m ' 


Ir s . 

rr 


(33) 


Then, applying v to 's m+l = 1 v n , , we have 

y'J*i ' Is n , + i " ‘ s ™ 


or 


(34) 


Also, applying v to n S m+1 = we have 


‘(Vi) 


v( n s mtl ) 


= He 


- “s 
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which, by using Equation (34), becomes 


”s 


m+ 1 


n S 


(35) 


The i, multiplying both sides of the Adams -Bashforth predictor and the Adams-Moulton cor- 
rector by v 1 , and similarly multiplying both sides of the Stormer predictor and Cowell corrector 
by V" 2 and using identities (34) and (35) we derive the following summed forms: 

Adams -Bashforth Predictor Summed Form 


V‘ l v v ' 

1 


y 


1 


h 



n 



2 





(36) 


where 


and 



j = 0 


where 


Adams-Moulton Corrector Summed Form 


= h | a o Is m + (% + O y» + ^ 


(37) 


and 


l 


* . y-* a i*-i 

a * i_> j + i 

i=o 
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Stormer Predictor Summed Form 


7" 2 V 2 


y 


m+ 1 






(38) 


where 


l 


and 



1 = 0 


Cowell Corrector Summed Form 


7 ~ 2 V 2 y 
J m 





(Po + £,*) ‘S m 


n 



i = 3 



(39) 


where 




i 


and 



j=0 


The meaning of l s and ri S m+1 can best be seen from their positions in an extended difference 
table (Figure 7). Examination of this table shows that the sums can be maintained by relationships 
(34) and (35) 


* S . , = 1 S - y ' , . 

m+1 m -'m+1 


' S m+ 1 = IIS m + * y^.+ l 
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I 





3 


^ rr» - 3 


V 2 y 


2 


1 m - 2 


V3> 


! 





Figure 7— Extended difference table showing I S m and n S m . 


but that initial values for some , s iB and 11 must be supplied. These initial values can be deter- 
mined by inverting the corrector formulas ( ! S m is eliminated from the Cowell corrector since its 
coefficient, p* + p*, is zero) and solving respectively for I S m _ 1 and II S ra ^ 1 


y'm-l 


- — v ' 1 + r + cL i*^ 2 y M i + • -1 

2 Vm-l a 2 3 *m-l I 


y rr — 1 


•n-« h [12 


n y -'-i + Vy ^'-i + + • • ■] 


(40) 


(41) 


D. Ordinate Forms 

All of the foregoing formulas involved difference operators. They are thus known as the 
difference forms and summed difference forms. Another useful form of these formulas which can 
be used under certain circumstances is the ordinate forms. 

When using the difference forms, the order of the method can be dynamically changed as the 
problem dictates. That is, on the basis of the number of corrector iterations, the order of the 
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method (the number of backpoints) could be increased (or perhaps decreased) to improve accuracy 
(or lower computation time). However, in satellite orbit determination, the functions are usually 
smooth enough so that the order of the method can be fixed. This permits us to take advantage of 
the ordinate forms of the Cowell and Adams type formulas. 

In using the difference forms, it is necessary to maintain a table of backpoints and tables of 
differences. The ordinate forms enable us to rely solely on the table of backpoints thus obviating 
the computation and maintenance of the difference tables. This simplifies the integration process 
and enhances calculation speed. 

Consider the Adams -Bashforth predictor ( 21 ) substituting definition ( 2 b) for 



Expanding the expression in brackets and denoting by z , we have 

a o(-D° (o) Z 0 + 

(o) Z o + c i( (i) z i 

(o) z o + (i) z > * M-i ) 2 (2) z * + 

V" 1 ) 0 (0) Z 0 : a 3 C" 1 ) 1 (1) *i + a 3 (-l) 2 (2) Z 3 + a 3 (- 1) 2 (3) z 3 - 


a 

n 


(- 1 )° 



a n (-l)’ 






0 
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Then collecting the coefficients of like ordinates, the expression becomes 




+ Z n-.(" 1 ) 


/ n- 1 

\ ( n 

\1 

a n-« \n - I 

1 +■ a [ . 


/ n \n - 1 

/] 


Zn(-l) 


" -(")] 


or 


y»M = vl + y; fe) - C, £}. (l) + y ;: 2 (i) 


y »-^ ZX (3) “ • + y ”-- 1 ' Z a ‘ (n- 1) + y -n (n) 


which can be represented as 


y 


m+ 1 
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where 


a 

) 


(- 1 )' 


£:■■(!) 


Sample calculations of the coefficients a, for a fifth order Adams-Bashforth predictor are given in 
Table 5. In like manner, the ordinate forms for any order of the summed and non-summed Cowell 
and Adams type formulas can be developed. 


Table 5 

Coefficients for Fixed, Fifth-Order, Ordinate Form Adams-Bashforth Predictor. 


-/o\ 

, /i\ 1 ^ /: 

!\ 5 , ( 

3\ 3 (< 

t\ 251 " 

1901 

\0/ 

1 \0/ 2 (c 

>/ 12 

0/ 8 it 

) j 720 

720 


17 M 

1 (■ 

l\ 5 ! 

3 \ 3 t (‘ 

4\ 25l" 

_ -1387 


2 (i 

1/12 i 

1 / 8 i: 

i 1 720 ^ 

360 



' / 2 \ 5 / 

3\ 3 

M —] 

^2 = 

_ \ 2 ) 12 + \ 

2/8 \: 

2 j 720 J 


, [(■ 

3\ 3 (‘ 

*\ 251 1 

C7 

3 = CD 3 |_i. 

3 ) 8 + (: 

3/ 720 J 


4 / 4 \ 251 ] _ 251 

a 4 = 4 / 720J “ 720 


Thus, the ordinate forms for the non-summed integration formulas are 

Adams-Bashforth Predictor Ordinate Form 


n 


1 = 0 


where 

- = (-o' £/■(!) 

i = j 


( 40 ) 


Adams -Moulton Corrector Ordinate Form 




+ h 


r-o 
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where 


C7 * 
J 



where 


Stormer Predictor Ordinate Form 


n 



i - J 


( 41 ) 


(42) 


Cowell Corrector Ordinate Form 


2 v„-i ■ y*-2 + h 


n 


j=0 


where 


n 



i = J 


(43) 


The coefficients , o', \ , > * are given in rational form in the appendix in Tables 5 through 
8. Within each table, subtables are presented on the basis of n ' 0 to n = 15. 

The summed ordinate forms are 

Adams-Bashforth Predictor Summed Ordinate Form 


y.n = h | a o ,S m + y " L i 


= (-»>* 2 ] a; (!) 
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where 


(44) 


Adams -Moulton Corrector Summed Ordinate Form 



1 (- 1 )' 


£■.-(!) 


where 


and 


a* - a * + 1 for i > 0 


( 45 ) 


Stormer Predictor Summed Ordinate Form 



n 



where 




( 46 ) 


Cowell Corrector Summed Ordinate Form 


y 


in 



V. 
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I 



(- 1 )' 


where 




i = J 



and 


(47) 


The coefficients >/, and \ * are given in rational form in the appendix in Tables 9 

through 12. Within each table, subtables are presented on the basis of n - 0 to n - 15. 


REMARKS 

In determining the orbits of artificial satellites, in which the equations that describe the 
satellite’s motion are extremely complex, numerical integration methods are very fruitful. 
Predictor -corrector methods for numerically integrating ordinary differential equations are used 
because they are efficient and lead to accurate results. In general, predictor -corrector methods 
have the following advantages: 

1. Generally only one or perhaps two evaluations of the function need be computed at each 
step of the integration whereas one-step methods require at least four or more evalua- 
tions of the function. 

2, The difference between predicted and corrected values provides a measure of the error 
being made at each step of the integration. Thus this error, which is better known as the 
local error, can be used to control the stepsize employed in the integration . 

Some disadvantages in using predictor -corrector methods are: 

1. The process is not self-starting. 

2. The process is highly complex to program. 

m 

The main sources of trouble that arise when using any type of numerical method for integra- 
ting ordinary differential equations are (Henrici): 

1. Truncation error due to finite approximations for the derivatives. 

2. Propagation errors (instability). 

3. Round-off errors due to a finite number of decimal figures used to express the coeffi- 
cients in the formulas. 
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APPENDIX 


The formulas for the coefficients presented in the following tables were programmed in 
fortran using a rational arithmetic package to eliminate the deterioration which would have been 
incurred using floating point arithmetic. This rational package consisted of the following subrou- 
tines: 

(1) GCD - A function which uses Euclid's algorithm to compute the Greatest Common 
Divisor of two numbers. 


[Aj , A 2 ] = GCD > 0 


where GCD = i if A 1 - o or a 2 = 0 or if a l or a 2 is not integral. 

(2) ADD - A subroutine which performs rational addition defined by 


Nt 



1 \[<VD 2 ] ) 2 V [D lf D 2 ] 


[ E r 


N 3 

[ D 3 .N 3 ] N 4 

d 3 = 


(3) SUB 


(4) MPY 


A subroutine which performs rational subtraction defined by 

N, N 2 N, (-Nj) N 3 
' D, = D, + ' D 2 = “ST ' 

A subroutine which performs rational multiplication defined by 

N| N a 

ll A [ N r°j ' [ N 3 D i] N 3 

Di D 2 D x D 2 D 3 

vw\ ’ 


(5) GRBC - A subroutine which calculates the Generalized Rational Binomial Coefficient 
defined by 
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s - (l - 1) 


where 




m - 0, 1 , 2. 


s r '2. -1. 0. 1. 2. 


and 


0 for m > 0, 



1 


(6) HS - A subroutine which rationally computes the coefficients of the Harmonic Series 
defined by 



1 = 1 


These subroutines were so constructed that the numerator and denominator of any result were 
relatively prime (i.e. (N, D) _ l). Also, the sign of any term was carried by the numerator while 
the denominator was kept positive. A zero denominator was used to indicate loss of integral sig- 
nificance in the computation of a term. 

These subroutines were used by a main routine to calculate the coefficients of the difference 
forms of the Cowell type formulas. A subroutine was used to calculate the coefficients for the 
ordinate forms. A final machine language subroutine was used to format and print the coefficients 
in rational form. 

Tables 1-4 give the coefficients of the difference formulas. The coefficients for the summed 
difference formulas are not presented since they can easily be taken from the non-summed 
coefficient tables. Tables 5-8 present the coefficients for the non-summed ordinate forms of the 
formulas. Tables 9-12 give the coefficients for the summed ordinate forms. Although the lower 
order ordinate forms are essentially meaningless, they are included in the tables to provide com- 
pleteness. 
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Table 1 


Table 2 


Adams -Bashforth Predictor, 
Non-Summed Difference Form 


Adams -Moulton Corrector, 
Nori-Summed Difference Form 


1 1/2 
x 2 / 1 2 
a 3 ?>/b 

X A 2 S I / / 2 J 

v s / 2 a a 

*6 1 9 (J« 7 / 6 csao 

1 7 6267/ 1 7 2 fi G 

*8 l 0 7 OU 1 ; / J 6 2 a 8 nr, 

*9 2->7 13/89600 

1 io 2*642263/96000*20 

1 11 H777223/17S107MO 

1 12 7^3604264367/26163487 *> aUDG 

*13 136364763817/40236134 *i n 0 D 

1 14 1 2 3 I * 6 J 0 9 8 1 9 6 6 7 /44034S4976DUQ 

1 15 2622 I 4^6/98402304 



a* 


a s 

a 9* 

a io 

a ll* 


a !5 


1 

- 1/2 
-1/12 
-1/24 
-19/720 
- 3 / 1 6 J 
-8*3/6 OHflG 
-276/24 l V 2 
-3396^/3626000 
-0183/10366 no 
- 3230433/47400160 n 
-h 67 l / 7 0 d 4 B 0 

-l 3^9677909 3/26163487 3 6GnO 
-2224234463/4736] 7962000 
-1 322B2G4J177/3 1 3 H 4 184 m 3 2 000 
-26396blO63/6097623OsOO(J 


Table 3 


Table 4 


Stdrmer Predictor, 
Non-Summed Difference Form 


Cowell Corrector, 
Non-Summed Difference Form 


0 

1 

2 

3 

4 

5 

6 
7 
3 
? 

10 

11 

12 

13 

14 

15 


1 

G 

l / L 2 
1/12 
1 V / 2 4 J 
3/40 

8 6 3 / J 2 U 9 6 

2 7 i / 4 0 3 2 

3 3 9 6 3/6 1 8 4 0 0 
6183/129600 
32604 33 /6 32 22400 
4671/74848 

1 3 o 9 % 7 7 9 0 9 3 / 2 3 7 7 6 R 9 7 * 0 0 0 
2224234463/39626496000 
I 3 2 2 6 2 8 4 0 1 2 7/24 14168064 U 0 0 
263 9 6S 1063/49268736001) 



l 

- 1/1 

1/12 

0 

-1/240 
-1/240 
-221/6U46G . 

-19/6048 
-9829/3628800 
-4U7/I72800 
-33J l 67/ 1 69667200 
-24377/1 3306600 
-4261 ) 6 4 4.7 7 / 2 6 I 5 3 4 B 7 3 A 0 1]0 
-7 UO 7 4463/ 4 766 1796200 
- I l 9 7 62206 7/89 6690996 200 
-97997961/80472268800 
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Table 5 


Order 

Order 

Order 

Order 

Order 

Order 

Order 

Order 


Adams -Bashforth Predictor, Non-Summed Ordinate Form 


2 

°0 

3/2 


a l 

-1/2 


3 

a o 

2 3/12 



-M/3 


CT 2 

5/12 


4 

°o 

55/2M 


CT 1 

-5^ / 2 M 


CT 2 

3 7 / 2 M 


CT 3 

-3/8 


5 

^0 

1701/720 


a i 

- i 3 o 7 / 3 6 0 


^2 

109/30 


^3 

, -637/360 


^4 

251/720 


6 

a o ■ 

M?77/lM'4u 



- 2 6 M l / H 8 Q 


a 2 

4791/720 


a 3 

-3699 / 7 20 


*4 

7 S 9 / M A 0 



-95/233 


7 

CT 0 

1 9 H 7 2 1 / 6 G '4 8 a 


CT 1 

- 1 ti A 3 7 / 2 5 2 u_ 


^2 

i JS 1 H~3/ 2 0 1 6 0 


a 3 

- 1 075H /9M5 


a 4 

1 3 5 713/20 160 


CT 5 

-3®03/252u 


a 6 

1 9 0 3 7/ 60 M 8 U 


8 

^0 

1 6 0 6 3/ MS 0 u 


CTj 

-l 152169/12096 0 


CT 2 

2M2653/ 1 3MmO 


°3 

-296053/ 1 3 H h 0 


a 4 

2 10 2 29 3/ 1 2096 0 


^5 

- 11 5 7 97/ 1 3MM0 



3 2 36 3/ MM MO 


CT 7 

-*257/17280 



Order = 9 


^0 

l M0<» /2M7/ J62BB00 

CT 1 

-215o2a03/ Id 14 4 0 0 


47*7 38 393/ 1 8 1 44 00 

a 3 

-69927631/1814400 


862303/22660 

CT s 

• 45586 3 2 f/ l 3 J 4400 

^6 

1 94 I 9743/ 1 a 1 4400 

a 7 

" -4032 353/ 18 1 440 Q 

^8 

1070017/3628000 


Order = 10 


4 32S32 l / l 0 16 8 0 0 
• i' 04995109/ 7 2 5 7 6 0 0_ 
66 40 3 1 7/ 1 H 1440 
- 2 8J4J 6 J 6 1 / 8 5 3 bOi) 

” 5 6 9 i if l 9 l 9 / 3 6 2 6 B 0 U 
-72234608 I / 3620800 

1 S 7 9 6 6 3 9/4 S 3600 
-23S7683/IM1440 

2 0 4 8 4 ri 1 I / 7 25 7600 

-25713/09600 


Order = 11 


2i32b09SA7/479GOl600 
-2n6794B781/l 1 97SUH0O 
! c 7 2 7 3 7 5 B 7 / 3 1 9 3 34 4 L 
■x io2 1 376209/ 1 9 9 5 & 4 0 0 
fr. 39 7 948 3 1 / 2 6 6 1 J 2 00 
-4226 u 679/623700 
? j 9 2 0 6 4 9 1 3/266 l 1 20u 
- i 0 6 0 8 0 2 9 l / 3 9 9 1 600 
“ 2 i 7 2 6 3 4 3 1 7 / I S 9 6 6 7 2 (, G 
-520419*1 /l 7 1G72UO 
268 4 2 2~S 3/95 0 0 0321' 


Order = 12 


6 

a 7 


4*27766399/958003200 
-6*779 3 4721/3 19JJ4400 
12126645437/191600640 
-1 Sn64372V7 3/ I 0 6 4 4 4 8 0 0 
35*4987256 1 /1S96672U0 
-4l7?U273229/|S9o6720u 
3 S 1 'i 3 9 2 0 4 8 3 / 1 S 9 6 6 / 2 Li 0 
-*>5551 749/436 1920 
92363662 9/1520640 0 
-1 7 4 LQ24027 1/950OO32UO 
30082309/9123840 
-477 7 223/1 741 8 2 40 


Order = 13 


] 3 U 04406523627/26 I 53*073600 0 
-911701 102989/34626496000 
S 9/, 1794194517/726 4 0576 u 00 
-10498491594103/52306974720 
207 .Li 767 690 l 3 I / 58 1 i 0360000 
-3476616/9 IS 0 49/726 4 857b 0 00 
7 2 41 3 3 0 14 5 3 3 / 4 8 6 4 8 6 U (j 3 
-2 02690 rj 57 7 631/0 0 7206Hu0Ij 
2 2519 5 7170793/1 1623772loU 
-20212291373037/261534873600 
4 5 n ;1 4 1 4 5 5 5 2 0 l / 2 I 7 9 4 S 7 * M □ 0 0 
-169^39834421/4 H 432384 u 00 
703*04254357/2615348736000 



Order = 14 


7 

J 8 
7 9 
7 10 


<J 

<J 

<J 


n 

12 

13 


o057j 0705/17220403? 

- ! 4 09 7 1 1 7 SO 6 7 9621/5230697 472000 
8 9 S 4 l 1 7 5 4 1 V 7 7 7 / 8 7 1 7 0 2 9 \ 2 0 3 U 

- 344 i 2*22659 J9 3/ 1 746HG4 *6000 

5 7 0 A d 5 Q 1 H 3 5 d 16 1 / l J'l b 1 J 9 4 9 4 M G 0 
-3 1457^3595 JH13/38745 9 072 Ofl 
i3HnM6u25652H57/MS297lb2U0U 
-350379^27127877/435091456000 
3 l 0 4 7 9 q 5 5 ft 7 5 4 5 3 / 5 0 1 186608000 

- 1 0 3 2 0 7 8 7 6 O h i 3 / 3 8 7 H 5 9 0 7 7 0 0 

7777659iS*9H9/7H72H2HV600 
- 210291621 13651 / 6717^2912000 
6440951 1 97929/1 7435*5624000 
“10 4 1647 6 3817/4 02061 3 4 H U 0 U 


Order = 15 a Q 



1 3 3 2 5 6 5 3 7 3 n 3 7 3 / 2 4 l 4 1 6 8 0 6 4 Oil U 
-600076741 50257/196151 15520 0(J 
39Aii4 2l*.7021j4fM/3|3iS4l84rtJ2£JU0 

- 7 59 9 1262345039/ 700^398 4000 
2529rt4l0i37uMl4?9/Jl3rj4164H320 0U 
- 2 6 1 4 0 7 9 .1 7 LI 7 8 i 7 3 3 / 1 9 6 i 5 1 15520UU 
170236 7 55 5 33] J 5 L 3 / t U 4 6 1 394944000 
- 2 1 * * a 15342637 / ) 2 7 70257 5 0 
1 3 7 6 0 O 7 7 I 1 2 0 9 4 7 5 3 / 1 u 4 to J 3 9 4 9 4 4 ?JUO 
- 154403 i 4 744 7 * 463/ 1 9 6 J 5 1 15S20JU 
1600835*79073597/4403454976000 
-532^2613384023/490377868000 
8592364766 0 -»231/313il 4184932 00 0 
-69 646 1 4 42637 / 1 783 1 92 32000 
1 l ft 6 3 0 9 ft I 9657/44834549 7 6 000 


Order - 16 



36255St26427073/b7768369664000 
- 2161567*7124 5 849/6276636 9 664000 
.7 4016 I 3aU731V49/4n2833612B0u0 
-4372441680074367/89869 0 99 5 2 QUO 
72S53U 7 572759733/62768369 66 400 IJ 
-1319631 9 194J623 5 81/6276 6 369664uUu 
624877 ; 3i7o987o3|/2M922789*60OOO 

- 7 t) n U 6 m 6 79 / 0 7739^3/20922789868 0 00 
6 2 0 2 9 i m l a 2 I 190 4 81/2 0 922789rtH8uOO 
- 1 299 30044 i a 4 2 J 7 33 1 / 62 7bfl 369 664000 
1 0 l U 3 4 / 8 7 ? 7 b 4 7 0 * 9 / 8 9 6 6 9 0 9 9 5 2 U 0 U 
-2A 7 4 358uj7386529/57062t54740Ufj 

903857 1. 7 52734 tJ87/b;>76H3bV6o40u 1.1 
-17344 4 3196692599/62760 3 69664UU0 
3 o 8 0 7 1 4 2 2 7 3 & 6 9 / ft 9 6 6 9 0 V V 5 2 0 0 U 
-26221445/98402304 
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Table 6 

Adams-Moulton Corrector, Non-Summed Ordinate Form 


Order = 1 a* i 


Order - 2 

CT 0 
rr * 

1/2 


a \ 

1/2 

Order = 3 

a o 

S / l 2 


a l' 

2/3 


a 2 

-L/12 


Order = 4 

a o 

3/8 


a r 

19/24 


^2* 

-5/24 


^3 

1/24 


Order = 5 

CT 0* 

25 l / 720 



323/3*0 



- n /3o 


a 3 

'-33/363 


°4 

-19/720 

Order = 6 

*0 

9 5 / 2 d 8 



t 4 2 7 / M 4 U 



- 1 33/240 



24 !•/ 720 



-173/L44u 


a i 

3/160 

Order = 7 

°0 

1 Vt)8 7 / 604BD 



2713/2520 


CT 2 

- 1 5 4 M 7 >20 16 0 


a 3* 

Se)6/945 


a 4* 

-V/37/20lo0 


a s 

263/2520 


*6* 

-.363/60480 

Order = 8 

CT 0 

5257/17280 


°T* 

1 39649/120960 



- 4 S 1 1/4480 


a 3* 

l 2 3 1 33/ l 2 0 9 6 Q 


a 4* 

- R 3 5 4 7/1 2 0 V 60 


a 5 

1537/4480 


CT 6 

- 1 l 3S l / 1 209611 



275/24 1 92 


41 



Order 


9 


Order = 10 


a o 

l 0 7 0 0 1 7/3625900 

a i 

2213547/1814400 

a 2 

-2302297/ 1 8 l 4HQ0 

a 3 

2797#, 79/18144 00 


-3 l 957/22600 


1573169/1014400 


-64*6 0 7/181 4 4 U 0 

a 7 

150437/1814400 

a i 

-33953/362880D 


°- 0 * 

2 5 713/89600 

CT r 

94497 1 7/7257600 


-1406913/907200 

^3* 

2 0 0 0 7 9/90720 


-9641823/3626600 

ff S* 

675504 1 /3626800 

a 6 

-462127/463600 


336963/907200 


- 1 l 6o87/ l 45 1 520 


8 183/ 1 036800 


Order = 31 cr 0 * 

<7* 



26B42253/9S80G32L 
*6404^413/1 t 9 7 5 0 4 u 0 
-796726163/159667200 
1206 1 7D9 / 39 V 1 6 3 0 
-33765029/8870400 
222/571 /62 3700 
-21677723/0870900 
236 4379 I / 1 99&840D 
-123169 13/31 9 3 3 9 9 U 
907 1 2 \ 9/ 1 i 9 7 S 0 9 (j fl 
-3 250433/479 001 600 


Order = 12 



9777223/179 1 82HU 
I I 79 79 92 1 9/9 S 8 (JO 320 0 
-996929 1 3/956 1 920u 
36465037/9 1 2 J fl 9 0 
-102212233/1/ 7 9 0 8 0 U 
in 07253 bhl /1596672DD 
-9191 0991/17790800 
Q 0 l 2 8 9 9 H 3 / 15 9 6 6 7 2 U 0 
-G7U6479 1 763966880 
>8970 93 27/9580032 U0 
-68928781/958003200 
4 6 7 ) / 7 8 8 9 8 O’ 


Order = 13 



7 0 3 6 0 9 2 5 9 3 b 7 / 2 6 15 3 4 073 6.0 0 0 

A c 95 2 0 90 6 9/490294 4 0 Q U 

-55 1 1684 1 3 1)9/217945728000 
1 3 4 6c 7 7 4 2 565 1 / 261 534873600 
-4 05c. 0084533 1 /58 1 1 08600 00 
8 4400rt3~5 409/807206 40 00 
-4B 74320027/486486000 
S29 3940459 1 1/72648576000 
-229682444333/581 18860800 
4 0 6 3 3 2 7 8 6 31 7/26 1534873600 
-30336027563/72648576000 
27 24 8 9 I 25T/39626 4 96 000 
-13A95779093/261534B736000 



Order = 14 





12 


104*64763817/402361344000 
7 4 1 i 9708747 1 / 4 7 5 5 I 7952000 
l 6*8 ? 35 9 "4*5 3 7 9/ 5 8 U 8 8608 00 
16964 «‘9 5066809/26153 '4 8736000 
“17 4 8 > 4 A 00 3 7 5 1 / I 4 9 44 8 4 4 9 2 0 0 
957Sc8096b507 /58 1 188608000 
“ 7 R 6 6 I I 1554491 / 4 3589 1 456000 
1 3 3 5 fi I 701 7153/87] 78291200 
- 5 7 9 7 q H 56 53 6 2 9/58 ] J 6'8 6 0 8 GO u 
5 1 2 4 n 5 1955567/ ] 0461 39494400 
- 4 5 9 rift 1 7 « 0 2 5 67 / 2 6 "! 5 3 4 8 7 3*6 0 U 0 
1636*420501/3773952000 
-69r9l417279/l0He>l394944UU 
2?24234443/4755l7952000 


Order = 15 


a 2 


CT 4* 

<9 




11 

r u 

7 13 
T H 


1|66i09hi 965 7/H*483M5HV760U0 
3173 18SH7U92 9/ 196151 1552000 
“102**5 lHR9S 6 217/313*418*4832000 
39 137Ul*47a249/*490377883J00 
“ 71353086250691/4483454976000 
48 «69«*761294 77/196151 1552000 
- 32123140027491 1/10461394944000 
18029005269/ 1 27 7 0 25750 
-236770944732449/10461394944000 
26159487787579/196151 1552000 
-1575 0493659/ 931/31384184832000 
137*55863 I 53/7 0053984000 
-141 10460969927/31384184832000 
124922452271/196151 1552000 
-l 3 7782 “40 127/ 3 13 3 4184632000 


Order = 16 



2522 1 445/98402304 
105 14505875 7073/62768369664000 
-2099728761 1259/57Q6215424000 
6 l 2 7 4 4c 4 1 06 5 3 37/62 769 369 66 4 000 
-18954ft 3 904 3 6*67/8966909952000 
2285l6ftt9«347 7 33/62768369664000 
-3l2745in7l99J5fll/627b836V664U00 
1 13931390961 7631/2092278 9ft 08000 
-998781^7 6758233/20922789888000 
67978 1959*4688 1 /2092278 9 880000 
-1096355735402331/62768369664000 
64486 1584 19 '.-169/8966909952 QUO 
-13751 5 71378 9 319/627 4 8369664000 
277 ) 9184284037/62768369644000 
-346748936/599/62768369664000 
2639451 0 53/989762304000 
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Table 7 


Stdrmer Predictor, Non-Summed Ordinate Form 


Order =1 k Q 

Order =2 k Q 


Order = 3 A. 0 


13/12 
- I /6 
l/l 2 


Order = 4 


K 


o 


>s 

\ 

k 


1 

2 
3 


7/6 
-5/1 2 
\ / 3 
-1/12 


Order =5 \ ?v9/;hq 

-11/15 

^2 97/120 

^3 -2/5 

\ I9/2M0 


Order = 6 


317/240 
- 1 33/ 1 2 J 
1 « 7 / 1 2 0 
-23/2U 
1 C 9 / 2 H U 
-3/40 


Order = 7 


*0 8Ml99/hOHBO 

-lb*H7/100riQ 
529 * 1/20160 
- 3 3 9 3 / iS 120 
\ 3U731/2O1A0 

-i0 7j/ujaa 

8*3/12096 


Order - 8 


k x 

k „ 


22081/18120 
-HSl 1/22HJ 

n o 9 3 3 / 1 uaag 

-300227/6048Q 
9857 /252G 
-39017/20160 
33 1 9 /*048 
-2/5/9032 



Order = 9 


'j^Vl » 1/J62 <i8U0 

-23.J229 7 /90720U 
5347667/907200 
-;a3u7wV/90720J 
61 'i 621/ 7 2576 
-5 uHJI 6V/90 7 2 0 0 
2 1 & ; S *• 7 / 9 0 ; 2 0 u 
- 3 3 / iM 7 / 9 0 7 2 J 0 
3 3 1 b 3 / 5 I 8 4 Ll a 


Order = 10 


1 153247/725760 
-1 *4 [J 0 v I 3/4S3ftJU 
i 4 0 9 7 M 3 / V I J ; 2 0 J 
-12642403/907200 
2 V a 5 o 3 3 7 / 1 8 i 4 4 0 0 
-24601|3/lAl*i4j 
6 7 7 3 1 S I / 90 7 2 0 1 
- 2 b9 >333/907200 
3 2 <» 5 4 1 / 5 1 3 4 0 J 
-3 1 A3 / 1 2960 TJ 


Order = 11 


763465639/ 1 59b6 7 2UM 
•29o7?S183/79A3JoJu 
1 7 ‘42930/63/ 1 596672uU 
- 4 2 *1*4 0 2 3 5 I / 1 9 9 5 a 4 Q o 
2/3 73 *11223/7 9833600 
- J ! 55556697 / 19 9 1 6 8 00 
lA3752 3*»63//vi)J36U(j 
•29^64973/285 1 200 
* 39 9 9 9043/ 1 59o6 /2u0 
-5379/223/79833600 
3250433/53 222400 


Order = 12 


1 9 m 9 4 6 G I / 1 1404800 
.99 6 42413/228 fl 9 6 0 u 
*4 0 'll 1623/2H512ULI 
- 4 9 !S & 9 1 s> 6 A 3 / ! 596*7200 
7/a^2USr-7/s7a29 0Q 
-4496090419/79833600 
oSrio2S 177/1 99SA4DG 
•237451 7 1 1 9/79833600 
in&0 39HH79/798J360lj 
-62782 7 0/ 1 / 1 59667200 
8*4 f> 7 l/l 10800 
-Ho 7 1/78098 


Order = 13 


4 A 2 t_ l 5 5 4 7 i 343/26 183 4 B 7 j a Q 1 1 0 
-4817634131 19/1039728640 0u 
7 a3Sa 239S4H93/4355914560QJ 
”57 1 *03503363/ 1 30767 4 3 6 0 0 
i 4 9 7 7 1 QB7 1 199/193729S3A00 
-lA5t«SS20S449/ 18162194000 
3l 4 *964546373/311 3 S1Q4Q0Q 
-13 68797279699 / IH162144QG0 
161456197531/3874590720 
-In 456898207 0 1/ 653 3 3718400 
I967fl57329773/435H914660Q0 
-0 1 742398949/ 108972364000 
l 3695779093/237758976000 



Order = 14 


\ 

K 

K 

K 

K 

k 

k 

k 

k 

k 

k 

k 

k 

k 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 
to 
11 
12 
13 


6 4 1 1 3 6 * 2 j 043/173*212*8000 
- 1 4873594*379/2 9QS943U4JQ 
97 3 4 n 1 7 1 237^7/43509 1 456000 
-3?0 76^rt;S99ijH7/A53MJ7lHMQDU 
2 U 4 i y 7 3 9 4 6 I 3*9/ l743565624Qu 
-2530/004074409/ 1 *5297 1 S 2 0 0 0 
0MJSAl533724L/4J56V|4b6OO 
•29 A 967J98357/ I 72 9 72300 0 
] 1 "i 3 3 pH 6S5 8 073/9686476*000 
- 1 4*3*921/ I 3 7 0 1 / 2 A l 5 3 H 9 7 J 6 0 0 
I 11*3270 3 09|/53 3137600 O 
-4*710*55301/3/1 78291 20 
M/trt6 0 6;2n7/237/S89/s00a 
-27242344*3/39626496000 


Order - 15 


'5 35773 9 661 13 3 / 2353137712000 
-10 24851 *395621 7/1S6V2 0 924160 (10 
3 4 3 2 2 ^9 331 12Ul/l2533673932au 
-125041 93071 17^1/15692392*1600 
54031 997UJ62267 1 /3lJb4104032uGO 
-4*5*63943*6 l 7463/1S6920924160JU 
3706/4 4 77952 <J 091/1 0 46J 394 9 440 UU 
-94 j 44^33621037/261 3 34373600 
582613*053861*7/20922789*88013 
-261 173190139471 1/18692 392*16000 
236687*1 22999363/31 30418*832UUU 
-I 9 6; 3 On 0 7 641 1 41 / 7846044205 0 00 
3*239*32144313/6276336566400 
•25 3 3707 0 5976] / 31 384 L 8 4 4 3 2 0 IJ 
1 372H244J 177/2* 14 16806 4 QUO 


Order = 16 


7877074749153/3923023104000 
-209972676I 1259/23 5 31077120 Ju. 
10 34 6 17^93 4 5973/3136418483200 
-14513674963751/139488265920 
1 9 7 5 8 4 7 3 7 2 6 l a 3 3 9 / 7 8 4 6 0 4 6 2 0 6 0 0 0 
-13958696412681209/31384184632000 
988 7964365**84539 / IS692G97* 1 6000 

81'4 9773547^S4] / I 3076743680U 
-136391 5 iA? 3 19 5227/3138 4 104832000 
3703l57fl2V222i23/l56920924l6ou0 
-3D82 10 9*27^3 3 329/3 1 384J 8*1032000 
15771040394 797/323069747200 
-4047002625 3 543/62760369664 00 
IO34?13182 l}41/12O7Q04u32OOO 
-7439*51053/49268736000 
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Table 8 


Cowell Corrector, Non-Summed Ordinate Form 


Order = 1 


i 

Order = 2 


0 



i 

Order = 3 

^0 

1/12 



5/6 



1/12 


Order = 4 


1/12 



S/6 


S* 

1/12 


S* 

U 


Order = 5 

k S 

1 9/240 


x r 

17/20 


s* 

7/120 



1/60 



-1/240 


Order = 6 


0/40 



209/240 


K i 

l / 6 0 



7/120 


X 4 * 

- 1/40 


\* 

l / 2 4 0 

Order = 7 

V 

963/1 2 0 9 6 


k i 

4 9 9 9 / 1 0080 


k 2 

-7*9/20160 


K 3 

l 997/ 15120 


k : 

- 1609/2 11 160 


K S 

263/1 0090 


K i 

-27 l / 6 0 H o 0 


Order = 8 

k 0 

275/4032 


k * 

1 39 3 1 / 1 S 1 20 


S * 

- 7099/20160 


k 3 

41 1/1360 



- 1 i 4 / 7/60460 


k i 

29/315 


\* 

- 5 17/20160 


K i 

1 9/6049 



Order 


9 


Order 


Order 


Order 


Order 


K 

K 

K 

K 

K 

K 

K 

K 

K 


o 

2 

3 

4 

5 

6 
7 # 
8 


3 3953/5 i840o 
4 2 h 7 5 9 / 453600 
“8 1 629/453600 
l I i 4 3 / 2 a 3 i 0 
- 27933/72576 
1 1 J5*3/4b360Q 
- 2 J 0 l 7/226800 
3627 /226HQU * 
-7829/3626800 


10 


^ 0 * 


K 7 


i \ 63 / i 2 9 fe Q □ 
694999//2S760 

- 2 8 j |,U / 9 J 7 2 0 0 
53e>iJr»3/V072lJG 

- 6 1 3 3 9 3/707200 
980713 / 1*14400 
-5431 1/181440 

99431/907200 
-2/11/M3HOU 
4(j7 / I 7 2 H 0 0 


11 



32SJ4 3 3/53222400 
312^0 27/3193344 
-57124721/15966/200 
167457m 1 / l 9 958400 
- 86645069/79333600 
4237 ^ 577/39916800 
-2342533/3143344 
713 78J7/199584U0 
-14674 153/ 159 66 7 2 On 
1 4338 i 9/79833600 
- 3 .1 i J 1 o 7 / 15966720 0 


12 


K s 

K' 


4671//8843 
79704557/ /98335U0 
-7321 7/41/15466/200 
45553097/39716800 
-l 3691 1^29/7 7 3 3 3 6 00 
761 62079/394 1680 u 
-1 261 33369/7983360J 
460*61 3/4989600 
. 6 6 t 4 J M 3 / I 5 7 o 6 / 2 (J 0 
'4 3rt 3 229/79 8 33600 
-354/921/15766/200 
24377/13305600 


13 


V 

V 

ij. 

K 4 

K 5 


l lA 9 S 77 rO 43/237 / 5 J 976 nD 0 
221A832 5 3067/217945724 iVJ-J 
-246Q77242177/435691456000 
I 9 6 7 4 l 13 >U i 9/ j 3Q76 7 436800 
-<M 921 4 6/ 453/193 /2453600 
1 1440 0317337/ 3 632428 6000 
-4628549 3 157/31 135104 00 0 
4 2 n 4 8 5 5 I 4 9 7/ 36324286 0 00 

-7 3019 73J/93/173/295 3600 
A 3 ?6i53 4 349/130767436800 
-56774^335/7/43589 1 456000 
46 Bn 454 7 3 7 / 2 1 7 9 45720 000 

-'OBI 1 6 4 4 / 7 /2615348736000 



Order = 14 


K 

K 

K 

K 

K 

K 

K 

K 


0 

* 

1 

2 

► 

3 

% 

A 

s 

* 

6 

* 

7 


'8 

S 

K* 


272423 < 4 * 3 / 3 9 6 2 6 4 4 6 0 □ 0 
*336249*961 /‘iJ3VHlbiOu0 
t / 2 1 /VH r ,7?A-J0j 

1 ,i 7 * 3 t» 4 6 j 7 / 9 ✓ 8 7 l 7 8 2V l 200 
-^.\snJ7^3's3n3//!(SlbiH67 16 0 0 

4 9-4* 12/343V7/46&&4766J03 
-204 |B 7 6 3 J S g 9 / 3 6 3 2 4 2 « ^ j 00 
I,vi^2is^^i7/2I 794*7/6000 
-12977*94477/4 I S I 3 4 f 2 - J U 

2 mn* 7 26*297/1 7 43565 8 2400 
— 3 ♦» 3 7 3 S 7 7 <» I l 3 / o 5 3 ft 3 / 1 o 9 0 Q j 

S9MAi|6JSl9/M3S^9lH r 3feL!0J 
^nt,iij6/ r it«7/‘03^l V_, 6 0 Q 0 
7 0 0 7 i-4 A 3 / 4 7 * & \ / V t> 2 GO 


Order - 15 


s 

s 

S’ 

X * 
4 


S 

S’ 

A* 


i 322*2*4 0127/2414160064000 
3 3 4 1 6 3 U 4 6 2 0 t / 3 2 0 2 4 6 7 u *-* 0 0 0 
-73209721 1 39 ij 79/3 1 3IH 109rt32uOO 
374 1 in | 6 7 MM 1 7 la692J92sUuu 
-3 3ft *Vi 0*123*321/62/63 3 696*400 
l22 0 0M9b29^2339/l 5 S6V2J < ?2M16JaU 
• 1 )3 7 A 3 ? 9 4 7 9 OS 5 7 / I 0 9 6 l 3 9 * 9 4 4 0 U 0 
J 2 2 5 4 6 6 0 322 3 37/133 7 6743* ft QUO 
-2 l 3 3 71J6S1 M3l/2VM9o99d40J 
2o4346l7&4441 /6276ft 3 646640 
• !i 9 ? 7 -4 n 3 7 H 7 1 4 0 9 / 3 \ MM lft4ft320O0 
4 4 4 •4 7Ml23ftH37/7B4*0 4 h 2 Q ft 0 0 0 
- -hr, / ri 3921 3359/ 313 4 4 18 4 332000 
s3Xjft39fcR2)/ M3fe4lH4rt32ClG 
-1 l?7622Jft7/rt9469j99520u 


Order = 16 


\r 

s* 

\o 

s; 

*13 

*«; 

*« 


7*34651053/49268736000 
4 21H156H0 7^31/ 3 92302310*400 i) 
-29^17337232 3 29/313841B4332U00 
9 2994.56983377/313*418433200 
-412397630 7*291/6276836966400 
0 9 * 9 5 * 4 l 5 4 4347 / 784604620*000 
-5 4442«i:i3tt64569/34<l?13l64ft0JG 
4000973424 337 3/b 230697472 000 
- 31328932761427/2092273948300 

1 .M 0 ft 1 3 0 ft 3 7 / 4 ft a 7 5 5 7 7 0 
- 1 7 4 O 4 .'1 2 6 7 344 1 39/3 1 38 4 1 84832003 
35ft V<|6 7 u 602*9/15692 'J 92416000 

-214 *3 677 * 014309/3138 4 184832000 

232305003033/ 1*69209241600 

-l 2304 U4S72/9/', 276*369664 0 0 
9749790 l 7*04/2268800 
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Table 9 


Adams -Bashforth Predictor, Summed Ordinate Form 


Order = 1 

a o 

1 / 2 

Order = 2 

a i 

1 1 / l 2 


a \ 

-5/ J 2 

Order = 3 

a o 

3 l f l 4 


a i 

-7/6 


a 2 

3/8 


Order = 4 

a o 

i L ? i / 723 


a i 

- : 7 7 / tj *. 1 


a 2 

J 4 ] / 2 H 0 


a i 

-23 1 / 72 J 

Order = 5 

a o 

2837/1440 


a l 

-2343//20 


a i 

1 7/S 


a 3 

- 1 2 a L / 7 <? o 


*4 • 

95/280 


Order = 6 

a o‘ 

1 3^24 1 /60980 


a i 

- 3 Q 9 □ 4 7 / 6 0 4 6 0 


a 2 

19*251/30240 


a 3 

-1 45477/30240 


a 4 

23J 7 7/12J96 


a 5 

-19OB7/604H0 


Order = 7 

a o' 

1 1603/4430 


a l 

- t as 44 1 / | S 1 2 0 


a 2 

13^Vd9/] 20 9 A U 


a 3 

-20 i, 17/1890 


a 4 

l 5 6 5 5 J /2 4 1 9 2 


a S 

-32 37 l / 1 5 i 2D 


a # 

5257/1 7280 


Order = 8 

a o 

l 04684*4 7 /362d0 JO 


a l 

-326So 739/ 36238 JO 


a 2 

693JUO3/4QJ203 


a 3 

- i b4 0 7 0 47/725 7 A0 


a V 

12186649/723760 


a 5 

-3 3S9933/4U320 J 


a 6 

i 22 7 72 7/S 1 9400 


a 7 

-i u 70111 7/3628800 



Order = 9 


32 0 85 ? \ / \ 0 3ft8UQ 
-409R777 1 / J628800 
102 l 9 84 | / 4 U J 2 0 0 
- i 35352 3 I 9/3628800 
l 67287/4536 
-9x39609/403200 
5 3 9 3 233/518400 
-94fJU29/36268G0 
25713/39600 


Order = 10 


a 2 


7 


1 4535079x7/47900 1 6nn 
-22U609S7 I 9/ | 59x67200 
235733 Cl 09/6652800 
-4070AHS9 I /997«20Q 
I 152537553/ 159667 20 
-ltfl6893U 49/2661 1200 
-5376023/ 156400 
- 753022557 / I 99S64GC 
I 494 8 4787/53222400 
•26842253/95800320 


Order = 11 


a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 


0 

« 

1 

2 

3 

4 

5 

6 
7 
B 
9 

iq 


\ *69763199/ 9 58003200 
-396x01 17 41/23950 0 80 0 
1 *95 t 54823/3S4 8 160LI 
- I 7 4 7363563/ 1 330560 0 
4 1 1430596 1/31933440 
-4 9S967 /385U 
2o8/o 83637 / 228096CH; 
-86x56259/ 1900000 
7 6794*319 / 5068800 
- \ 4 4 7 9 4 7 5 9 / 4 7 9 0 0 1 6 0 
4 7 7 7 2 2 3 / 1 7 4 l 8 2 4 u 


Order = 12 


8 

a 9 

a io 

a ll 


I 0 4 4 9 n 5 7 7 3 / 6 2 7 / 2 6 I b 3 4 0 7 3 6 rj u U 
-5 104/4 4 950096 4 7/26153 4 8736000 
3 27?'* A 19198593/523069747200 
-2 40 68 n 96927 4 79/1743565 8 2H0U 
l 9 fl 4 6 0 2 i) 7 1 4 5 7 / 4 7 1 7 5 2 9 i 2 0 0 
• 1576 l abe 733287/622702080 00 
l 3 4 3 9 t 6 9 1 2 6 V 3 7 / 6 2 2 7 0 2 0 6 0 U 0 
-1 171404946 07U3/ 6 71 7 8291200 
1036 1259060489/ 174356582400 

-9 32na055662i)7/5230 69?H720U 
7 6Cal3348 307/2377589 7 60 DO 
-7 fl 3 /■ 0-25 4 357 /26 15348 7 36 DUO 


Order = 13 


7 

a 8 

a 9 


7335 2 £» 1 73 / 1 72204032 
-49;,M’J<j46587j73/26l534b736DUU 
104 6 39? 98352 2 9/i 3076743680 0 0 
-1026756 I 92 3 4099/5 23069 74 7200 
12lM4«n9lVo332i/34a7 131648GQ 
-4 0 3 IS 732897599/87 17829 1200 
3! 974 145463/69498000 
- 1 4 9 «, 3 1 3 1 9 6 5 5 0 l / 4 3 5 h 9 1 4 5 6 U Oil 
<>6 3 93nul79h» 4 71/34871 3 164600 
-1> 247 r, 42*72623/17435*56240 0 
4917039 137 l 7/23775897x00 
- 9 o P ri n 5 5 9 3 2 0 6 3 / 2 6 15348) 36 000 
lllol* 4/63817/402361344000 



'1091 m8*67*O73/2HiS16eQAM0U0 

-3ia?7jr,b263 7 263/il3hHl8**a3 2UD0 
S2^^9l3 529297u3/ , 3 23Q6V7n/2UUU 
-4?477ss70649V627/lb6920V2416Gu0 
A72|3a«3*9632H7/12S&3673932B0 
- 7 1 r\0 ? 0 ^ 'J H S 3 ft 0 6 1 7 / J 4 a 7 1 3 1 ft 4 8 0 fj 0 

3J8*KO*6*327bb9/373fc212M&GUu 

-20**s61g|7H2/6*3/2Mb3S8736UU0 
1 K 3 1 h n * i h 7 S M 3 A 7 / 3 s g 7 1 3 1 6 4 6 Out) 
■32^F-73*?331 3*7 l 7 / 1 2bS3673V32«G 
t i b ft 3 4» 4 2 ft 4 l 1 IjC 7/ l s26SS38b*U00 
“ > 2 4 i 3 4 7Jb2b29b3/S23u6V74720u0 
ft A ( ] 7 1 b 7 3 2 0 b u 1 / 2 4 IS iftBOAHOOLi 
-1 lg>>lCf9^l96b7/Ssa3H&SV7ftOOO 


29V7ft*7bA763G73/62768369fe6SUU0 
-1 1 o.U! >37 1553ft 1/3923023104000 
2bft*77l99&3MJIft7/2U9227B9ls88UU0 
-3 5<>VAb?79] 40 297/980/55776000 
jyn 8 ft 27 0/49/2510734/865*0 

-S71 1 9Si662 f JI»7s9/s3Sfc»9 l S b ft J 0 0 
SO 1 Ut.M 7 o 70s ? i ')9 7 / 2 9 88 9 ft 9 9 4 '4QU0 
-?t 3 2ibA3951 31/12770257*0 
90304MS7H7790HbV/697s26329 e 0uU 
-171s3iJl2«sit i 2o9/lbb92U92<f1b0 
2 Ci 0*5*5 4 7352 :j3b3/S7G62l542«*aUU 
“34 7 ftn 7 32493 03/297 i 9872000 
l 3 0 2 2 1 a 1 9 2 46*27/4828336128000 
-2lSf>«3S|\Rl l29/b6Q431B7 20JCi 
2bi 2 1445/98902304 



Table 10 


Adams -Moulton Corrector, Summed Ordinate Form 


Order = 1 

a o' 

1 / 2 

Order = 2 

a o' 

S / 1 2 


a f 

1/12 

Order = 3 

a o*i 

3/6 


a r. 

1 /6 


a 2 

-1/2*4 


Order = 4 

rr * ' 

25 1 / 720 


a *. 

S9/29Q 


°2*. 

-29/2*40 



1 9/ 720 

Order = 5 

a o*; 

95/260 


a i. 

7 7 / 2 *4 0 


J. 

•//Ju 


a 3, 

7 3 / 7 2 0 


a 4* 

-3/ 1 4*0 

Order = 6 

-o' 

19087/40980 


a r' 

2 3 / l 9 /60980 


a 2 

- 1 t 3 7 1 / 30 29Q 


- 

/3 i 1 /3D2iO 


a 4. 

-59*4 9 /60980 


a s* 

063/60900 


Order = 7 

a °. 

5257 / 1 7 2 6 Q 


a l 

6961/15120 


a l 

-66109/120960 


“3. 

3 3/70 


tt 4. 

“-3 15 23/ 1 2 09 60 


“5. 

t 297/ 15120 


a 6 

- 2 75/29 1 9 2 


Order = 8 

*0. 

1 G 7 JO 1 7/3628300 


a i. 

1 9093 1 1 /362d80D 


a i 

-299587/903200 


% 

1 15963/195152 



-926809/725760 


a s. 

1 1 29 77 /903200 



-273921/3628000 


-7* 

33953/3*20000 



Order = 9 


Order = 10 


Order 


11 


a o, 

a ? 

a ? 1 
a *' 

a 'c/ 

a,;' 


2b? I 3/87600 
m2 7967/725/60 
-39932 1 7 /3626B00 
S 30 52 / /H032DU 

-OH67/967U 
26 I 6 1 6 I /36 28*00 
-79019/80690 


a.* 

26 3 fi 7 7/3628800 

a l' 

1 

-d 183/1036800 

a c* 

26892203/95800320 

4* 

l03?9o939/lS96672uO 

a 2* 

-29 1 tb'393/l99S89QU 

a }[ 
a 4 * 

1 8H7 1351 /9979200 

-1 6931 9953/ 79833600 

a ** 

J 

a t[ 

a.** 

251 62927/15966720 
- 8 6 6 J 6 0 9 / 9 9 792UG 

63225/3/ 1 99S89QL 

v‘ 

s 

.1 101 1991/159667200 

325*19 33/979001 60 0 

CL » * 

0 

9777223/17918290 

a *’ 

l 

8 u 9 V H L> 1 / 1 190980 (i 

a * 1 

-6728320 9/956 1 9200 

a%‘ 

3 

1938. i 297/5702900 

< 

-07 2 631 0 1/159667200 

a r * 

76561 /299h8 

a 6 * 

-\372090l9/lS96672u0 

a„ # ' 

91021971/3991 6600 

< 

-107151937/31 93399UO 

< 

l 5 8 1 3 3 79/2 3 95000 HO 

a io 

-96/1/788930 


Order = 12 a () *' 



^□'5^n42S4JS7/26lblM87 360UU 
2QOt;ftQ673b3M3/26lb3H873A00u 
• 9 2 7 | 22699917/523069797200 

1 18068800959/39871316980 
-01079266099/67176291200 

j 9 1 7 9 9 H 2 9 □ 2 3 / 6 2 2 / 0 2 0 6 0 0 0 
-7*216313/0 3/622/070 Hu 0 0 . 

2 *4 0 7^^862 6 H 7/871 7 82912 00 
-8166391 l 0 5/ 6979263/96 

l8Slri999V/S9/b230697M7200 
“1661 97 0 939 73/261 5 3 4 8 7 36000 
13698779093/261 b 398736 OUU 


Order = 13 



1 0636 9 7638 l 7/9 0236 l 3 9 4 0 0 G 
*3 078 15 1 72993/373621296000 
-2 7 09 n 056 660 77 / 1 3076/ H36 8 000 
23 09 79 0 796 9 31 / 5 2 3 06 9 7 9 7 2 00 
-SO 799 2B3b99 3/69792632960 
9UP7n93u02299/9J58919S600U 
-22 1 SS33/2SQ2S0 
2«l 601 6S33S73/935 691 986000 
-I 7s 1 U 2 0 2 3 6 i 7/998 | 61 669U0 
1 9 9 6 9 0 9 9 :i 9 6 1 / 1 0 9 6 1 3 9 9 9 9 9 U 
-98*7720 70 77 1/1 307679368000 
I 6 099 8 2 8 36 S l / 2 6 1 S 39 8 7 3600 J 
-2729239963/97551 7982000 



Order 


14 


% 11663098(9657/4483454476000 

a i, 2504*31949133/2853107712000 

a 2 ( -12555699585959/5230697472000 

a 3* 83|9S>4054oij VI / 15692092416000 

a 4*‘ -6 46 3 130 I 3 3253 1/6 2 76036966 400 

a 5 ’ 50972790156553/3487 131648000 

a 6 t -4 2G7n«&74iim3/26l63487360Gu 

a 7*‘ 51 I6n779jb6'i 7 /3 7362.1 2460QQ 

a g ' -31 ! 73S.8779 135 1 /3487 1 3 1640OUU 

a 9 ' 27b9 79U2ti9S82|/6276836966MU0 

a io ( -247 S 77 l 1059413/1 5 6920924 16 IJUU 

*lt\ 20 4 0667423953/523069747200 0 

^12 * - ] 3 664 76 .3 96209/3 1 3 H 4 18 4 8 3 20 0 0 

132282840127/31364184832000 


Order = 15 



2522 1 445/98402304 
3 6 5 4 n 5 1 1 4 5 153/39230231 04000 
-I / 252? 3 H&4Q1S ij]/62768169664U0U 
2297797094;)8 3/32691B592OOU 
-88676146739 4 1 13/627*8369664000 
j49Anl7827JH9/l5692U924l6U 
- 1 923 3 84376 7 3 1 0 9/ 6974263296000 
3495 8 1 09 7 /1 J030B 75 
*4274ftQ96G6 1 69/9/20922789869000 
52S6nB2696499/435u9 1 466000 
-l 3579i 7 1932259/2510734786560 
1 743009 541 047/960755776000 
-8S1OA34387 4 37/2O92278V880GL1J 
2267l7a7 ill 1/39230231040GO 
-263965 MJ 53/609 762304L0U 
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Table 11 


Order = 
Order = 

Order = 
Order = 

Order = 
Order = 

Order = 7 
Order = 8 


Stormer Predictor, Sumined Ordinate Form 


! / 1 2 


K 

K 


0 

1 


l /a 
1/12 


**0 5 9 / 2 4 0 

^1 -29/120 

19/243 


7 / / 2 HQ 
-7/ 1 b 
7 3/240 
-3/40 


N ( 2 3 / 1 9 / 6 0 4 3 0 

*2 -1 13 7 1/15 12 0 

*3 / 3 8 1 / 10 0 « 0 

^4 -aSH9/ IS l^U 

S' 883/12096 




K 

K 

K 


3 

4 

5 


696 l / 1 S l 2U 
-6* 109/60480 
79/70 

-3 l 523/3CI24Q 
I 2 4 7 / 3 □ 2 4 
-275/4032 




1 904 3 1 l/ih28MiJ0 
-29 *58 7/201 6 0 0 
1 18463/48344 
-42*809/ 18 1 i40 
1 1 2 4 7 7 / a 0 6 4 6 
-278421/604800 
33983/818400 


^0 4 2 / 487 / 7 25760* 

* 1 - 3 4 V 8 2 I 7 / 1 8 I n 4 0 0 

*2 500327/134400 

S' -12939/2838 

S' 28 I 6 1 6 1 / 7 2 5 7 & 0 

S -24019/13440 

S 2630/7/5 1 9 400 

V -4133/129800 



Order = 9 


lti379J'U 9/1596672 u 0 
-241 IS *43/9979200 
14071351/3326400 
-157314453/19950400 
25162927/3193344 
-n 6 6 0609/ U»d32 JO 
*3225 73/2H5 1 2U0 
-1 t fJ I I4P1/L995840U 
325 JM 33 / 5 322240b 


Order = 10 


:V.J89hii 1 / l \ HiJ H 30j 
<•6 7 2 * 3209/22609600 
1'U 9 J 2 9 7 / l 9 0 0 3 0 J 
-O 7263181/19916600 
382605/24946 
-13 7 204 019/2661120 J 
M 1 J 2 1 M 7 1 / i 7 0 2 4 J Q 
- 1 J 7 1 S 1937/39916300 
1 So l J 3 ‘7 9 / 2 6 6 l I 2 0 U 
-4671/76648 


Order = 11 


\ 

K 


0 

1 


K 

K 

K 

K 

\ 

K 

K 

K 

K 


2 

3 

4 

5 

6 
7 
8* 

9 

10 


2 O 05 A 0 6 7 3 S 3 4 3 / 2 6 I b 3 4 6 7 3 6 0 U 0 
-9221 2264441 7/26153487*600 
l 1 * 0 6 ri 6 0 J M i 9 / 1 162J77216J 
-43 3H 7 926 6 8 4 9/ 2 1 794572400 
34 1749426 023/12HS40416J0 
-242J631374. 13/10378365 J 00 
/ 4 ]? 4 4 4 A 2 a 4 7 / 1 2 4 5 4 U 4 1 6 0 0 
- 4 * 6 6 4 4 1 l i j 5 / J 7 1 7 0 2 9 l 2 
l 4a l 6990 775 9/58 1 188608 UQ 
-in 6 14/ 043473/261534873500 
l 3 6 9 5 7 7 9 0 9 3 / 2 3 / 7 5 0 9 7 t> 0 0 u 


Order = 12 


o 

*3 

*4 

*5 

*6 

V 


K 


9 


K 

K 

K 


id 

u 

12 


30/5101 723M 3/37362 l 248000 
-2 7 0 9 fj 0 5 6 6 6 0 7 7 / 6 S 3 8 3 7 1 640UU 
23 0 97967469. 5 1/174356562400 
-SO 7642*35443/1 7435650240 
40 0/ 043002299/47178291200 
-6646599/ I 25 1 25 
2H i An 16533573/52270206000 
-175102023617/6227020800 
14 469094 59 61/1 1623772160 
-44 A7? 2 0 767 7 1 / l 30767436000 
|6 04 952536 51/23 7 7509/6000 
-2224234463/39626496000 



Order = 13 


2904631949133/2853107/12000 
-]255S 69v r >fib<JS9/26|53487 3 6GGU 
39)9 5MHS4 6u91/523Lj6974720iJD 
-A 9 631.VJl332S31/lS6V20924lAijU 
5 0? 7279 0 MASS 3/6 9 742A32960O 
-42tJ7.3qS/4bi Jl 3 /h3SA9m*60QJ 
Si |An7 7V,^6-37/S3 37S464uOO 
-JJt7.Jtri7791 J5i/-4 3S3 9t4b6UQJ 
27S9/9U239 ‘j <3 21/697 •» 2*32960 Q 
-24/5/71 1 OS 9 H| 3/1^69209241600 
20 9 366/H2S9S3/H7Sbl 7 9 S 2 0 0 J 
- I A 6 647639 A 209/26 t b J H A 73600 U 

[37?a2^HJI27/2HlHlodu6HJuu 


Order = 14 


A. 

\ 

\ 

K 

K 


A 

S 


\ 

K 

K 

K 

K 

\ 

K 


6 

7 

8 
9 
10 * 
11 
12 


K 


13 


3 A S H n 5 i Mb is 3/39230231 OHQOO 
-I 724? / 1 4 S 4 Q MO 1 / 3 1 3b 9 1 B 4 8 3 2 u 0 O 

2 2 9 2 7 9 7 w 9 4 u A 3 / 1 J 89729*4000 
•R86 7 6 | 46 7394 1 33 / l 56920924 I 6U 00 

3 n 76 rU 7 A 2 7 J i 9 / 3 ] 3d4|64832 

• i ? 2 3 3 4 4 3 7 * 9 3 1 09/1 16237 7 2 1 6 Q 1 1 il 
2 4 47067679/13 0 3 U d 7 a 
- 4 ? 7 4 3 9 qB{] A 1*979/26 I 5348736000 
S 2 M n 8 2 * 9 „ 4 99 / 4A4J2J84UG0 
-I 3 5 7 ? | 7 ]9322S9/2slJ7J47B6bA 
174**0 9 3 4 |iV( 7/3913 9616000 
-AS 3363 (38 7 437 / | 74Jb6Sb24Q0U 
7 > 6 7 1 7 b 7 U 1 I J / J rj 1 7 7 l u u A U 0 0 
-2639*51 'JS3/4926d73A0U0 


58 


i 



Table 12 


Cowell Corrector, Summed Ordinate Form 


Order - 1 

s*' 

1/12 

Order = 2 

s*; 

1/12 


s* 

D 

Order = 3 

V’ 

l 9 / 2 *4 Q 



1/120 


S*’ 

- 1 / 2 '4 0 


Order = 4 

S', 

1/4 0 


S* 

\ / id 


s*' 

-1/60 


S" 

1/2*40 


Order = 5 

S*‘ 

161/12096 


S* 

67/1390 


S' 

- A n 9 / 1 iJ 0 ii n 


S*. 

S' 

7 l / 3 7 8 u 


-221/60930 


Order = 6 

S*‘ 

2/5/4032 


s*; 

2 2 1 / 4 3 2 J 


S' 

-21 17/30240 


S 

253/5340 


S' 

- 1 17 1/60430 


S' 

19/6048 


Order = 7 

S' 

33953/5 1 8400 


S' 

40769/604500 


S 

-8 353 / 4 S 3h 4 


S* 

l J937/ 1 8 1 440 


S' 

- 1 4 5 1 3/24 l 92u 


s 

1 1 7 29/604 BOO 


S' 

• V d 2 9 / 362 ‘5 BOO 

Order = 8 

S*‘ 

* 1 A 3 / 1 2 9 6 U 0 


s* 

33 3 27/403 200 


s 

-9 <, 4 27/604600 


s* 

1 3 587 7 / 7 25 7 60 


S' 

-4307/30240 


S' 

53257/1 209600 


s* 

-14 '3 29 / 18 1 4400 


S' 

407 /17 2600 



Order = 9 


Order = 10 


*>■ 

3250*33/53222*00 

K i‘ 

5727*1/5702*00 

5* 

- 8 /01A A 1/39916800 

5*. 

* 0 2 b Ji 1/13305*0 0 

V 

-91 7 J 3 9 / 3 1 9 J 3 3 * 

5*. 

737 0.6A9/J99 1 A30u 

5* 

-1 02i / /9/1 33JSbQu 

5*. 

75*331 /3??iA80b 

> * 

'8 

-33 1 157/159 672^0 

i • • 

,v o 

3*7 1/788*8 

5 ;; 

2 l 2 l 5 3/ 1 H 1 **J0 

r ' 2 

-l I 33* 197/39716001; 

X * 1 

'"3 

6073979/13305600 

5* 

-H 1 3s ■•♦9 8 7/ 7 7 833600 

A * 1 
5 

663*07/ ! 5 9 o A 72 

'-s' 

-3 0 73**7/1 3 3(J5 6 0u 

5" 

33 3 7 0* 7 /3 9936800 

v 

-2962873/159*67200 


7-M77 / i 1 1 




Order = 11 



'V 


1 u9S77'JU'0/^U7bii97^ iGJ 
m 6 | 74^n9/; o li348'M6Q0 

i * ? 3 0 3 * 2 i) / 9 / / ) 79*572*00 
-I -1732-1*2**9 / ] ? 4 b *0*1600 
*092 0 -* 2 3 7 / •* 9*2080 U *1 

- 7 | i> 7 9 l fl * * 9 / I 2 * 5 * 0 * l 6 J 0 
*1 3 0 * 9 2 l .19/21 7 V * 5 7 2 1\ 0 G 

- t 76* b S') 0*3 / 19 3/2 

* 7 A ij J 1 8789/26153**73600 
-h? 8] l A * * 7 7 / 2 6 1 5 3 *8 7 3 6 Out) 


Order = 12 


/;?2*23**A3/39A26*96U03 

6ntS57M57/*C)2JAl3**ua 

- ,7 U 8 1 *333*7/65 3 b37 1 GHQ0 

?*799A*H69/2767S&*8UU 

- 5 8 7 A 1 l2J629/*J389l*5AOO 
9 195*9 0 997 7/8*22/9 / U A 00 0 
- 1 1 Q 3 0 7 A 7 * 1 / l 0 3 7 * 3 * n 0 0 

6 6 9 1 7 u 1 H o 7 1 / M 7 l 78291230 

- ? 7 7 3 7 0 0 O * M / a 7 1 7 8 2 9 i z 0 0 

1 7A 0522 9 231/! 7 * J 5 A 5 8 2 H 0 1) 
- i * a 8 J 9 7 7 / .1 3 3 v 7 s» 0 0 U 
7 ij 0 7 * * /, 3 / * 7 5 3 i 7 9 * 2 J 0 


Order = 13 


s* 


; "5 


l32?326*'J12//2*l*lA«iUA*00U 
» 9 2 * 1 3 0 I 7 / l | 6 2 3 7 7 / t A 
-18 17 06612*97/3*87 1 3 1&HH0J 
1 3 A .l 7 3 i 8 t S A 7 / | 1 2 0 8 o 3 7 * H 0 d 
- * * 7 n 8 9 a 9 .1 8 A 3 / 2 3 2 * 7 b I * J 2 J U 
7 A * 3 / 3 5 ■» 12 7/ ( J 37836* QUO 

i 3 U * * A i 4 79/790 5 9 * 30*0 
-77?20 0 5 (>377/7 7*9 181**00 
3 0** 1 57832*7/ / 8 * 6 0 * 6 4 (l e U d 
- * Q 0 J b H 3233/830269**000 
1 S 2 5 A 9 :, A 1 7/3/1 78291200 
-1 l9 7o22 0*7/3 9ot,9 0 99b200 



Order - 14 


7A39AS 1 i]S3/H9 2gB 7 3MIUO 
iaV733369 0l9/t(JMolJ9N?HHiJU 

<4 0n9/;?3 2 

^«?75lrtJl3rt9/3lJ.iHl8irt32U 
-6i-:?'i^B37S.*/A7/2.)92278 r /888uii 
2 M C* } Q 8 2 7 8 !> A * 7 / S 8 l l d 8 o <j <4 J Q - ) 
-239S13 0J ?'M! 1/S23 u69/N720j 
696SsSJ27i/J77v|HBH0J 
M 165SJ 7 99/2 3 2S7S 9832 00 
3V*Al3S99V/«3/JI3HHl68rU2U0 
-237^7168*41 7^9/ l i73»j6V7 1720JJ 
2|797?72l7/iV3>29!iJ6u0 
-in77SJ27fe9 73/A?7€)e36V66‘4au 
97997 9& |/ .4DH722AbeOD 


61 




